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Abstract 



■ The electromagnetic inclusive transverse response of nuclear matter at sat- 

in. 

uration density is studied within the correlated basis function perturbation 
theory for momentum transfers q from 300 to 550 MeV/c. The correlation 
. operator includes a Jastrow component, accounting for the short range repul- 

. sion, as well as longer range spin, tensor and isospin ones. Up to correlated 

ov 

. one particle-one hole intermediate states are considered. The spreading due 

to the decay of particle (hole) states into two particle-one hole (two hole- 
one particle) states is considered via a realistic optical potential model. The 

6 . 

Schiavilla-Pandharipande-Riska model for the two-body electromagnetic cur- 
rents, constructed so as to satisfy the continuity equation with realistic v 14 
potentials, is adopted. Currents due to intermediate A-isobar excitations, are 
also included. The global contribution of the two-body currents turns out to 
be positive and provides an enhancement of the one-body transverse response 
ranging from ~ 20% for the lower momenta to ~ 10% for the higher ones. 
This finding is in agreement with the Green's Function Monte Carlo studies of 
the transverse Euclidean response in A = 3, 4 nuclei and contradicts previous 
results obtained within the Fermi gas and shell models. The tensor-isospin 
component of the correlation is found to be the leading responsible for such 
a behavior. The nuclear matter response is compared to recent experimental 
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data on 40 Ca and 56 Fe. 
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I. INTRODUCTION 

The cross section for inclusive electron scattering at intermediate three-momentum trans- 
fers (q < 600 MeV/c) has been object of many theoretical and experimental investigations. 
In the one-photon exchange approximation the differential cross section is given by 

<il(1 a M { q ^R L { qi u) + [tan\ 9 -)-$-}R T { q) u)) , (1) 
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where Om is the Mott cross section, q 2 = u 2 — q 2 is the squared four-momentum transfer, 
9 is the scattering angle and Rl(t){q,^) is the longitudinal (transverse) separated response. 

The total cross section is indeed well described by a simple Fermi gas (FG) model IJ, but 
the agreement disappears when the longitudinal-transverse (L/T) separation 0-|J is carried 
out for medium-heavy nuclei. Rl is largely overestimated by the FG model. However, 
the quenching of the longitudinal response is now well understood in terms of short range 
dynamical correlations, induced by the strong nucleon-nucleon (NN) interaction, and of 
nucleon degrees of freedom alone. 

In Ref. |TIJ a realistic model of correlated nuclear matter (NM) was used to study Rl{q, 
at the NM empirical saturation density, in the framework of the correlated basis function 
(CBF) theory. The density dependent nuclear matter CBF results have then been used 
in Ref. f| to estimate the longitudinal response in 12 C, 40 Ca and 56 Fe in local density 
approximation (LDA). The overall agreement with the experimental data was shown to be 
satisfactory. 

The present understanding of RT{q,uj) is more uncertain. Recent experimental L/T 
separations in 40 Ca |6||| have provided a transverse response lower than previous estimates 
JD at the quasielastic (QE) peak. On the other side, theoretical realistic calculations (feasible 
in the longitudinal case) have been so far prevented by the complicated structure of the 
transition operator, containing one- and two-body currents. In light nuclei (A=3,4) only 
the Euclidean transverse response, with the full current operator, has been computed JTTJ 



using the exact Green's Function Monte Carlo (GFMC) technique and the realistic Argonne 
t>i4 NN interaction [12]. GFMC cannot be presently adopted in heavier nuclei, so studies of 



RT(q,uj) in these systems either have considered only the easier to address one-body piece 
fj~3| , [I"4]| , or have treated the two-body meson exchange currents (MEC) within independent 
particle models (IPM) fl5|--|T7f . The MEC were found to substantially increase the one-body 



response in Ref. JO], whereas the IPM calculations of Ref. [[17] point to a slight reduction. It 
is worth noticing that the aforementioned latest heavy nuclei separated Rx(q, oS) show a good 
agreement with theoretical responses containing only one-body currents §| , downplaying the 
role of the MEC and in contrast with both the light nuclei case and the old 40 Ca data. 

Aim of this work is to use CBF theory to compute the symmetric nuclear matter trans- 
verse response in order to ascertain how it is affected by the NN correlations. Particular 
attention will be devoted to their influence on the MEC contribution. The results presented 
in the paper have been obtained within the exchange current operator model developed by 



Schiavilla, Pandharipande and Riska (SPR) in Ref. [0. The SPR model satisfies the conti- 
nuity equation linking the current to the NN interaction and contains intermediate A-isobar 
excitations. For the sake of comparison, also the standard one-pion exchange currents have 
been used. 

In nuclear matter, CBF calculations are based upon a set of correlated wave functions 



\n) =S 



i<j 



\n) FG , (2) 



obtained by applying a symmetrized product of two-body correlation operators, f(i,j), 
to the FG states \n)pG- An effective structure for f(i,j) is 

f(iJ) = E fMfaW® , ojfW = (1, a t ■ aj, <g> (1, n ■ Tj ) , (3) 

9=1,6 

Sij being the tensor operator; this operatorial dependence resembles the structure of 
the NN interaction. Additional spin-orbit correlations are often considered in realistic NM 
groud state studies, however, they will not be taken into account in this work. f(i, j) depends 
upon a set of variational parameters, which are fixed by minimizing the expectation value 
of a realistic, non relativistic hamiltonian on the correlated groud state. The g.s. energy 
is calculated via the Fermi hypernetted chain (FHNC) cluster summation technique JT9 



and the operatorial contributions (q > 1) are implemented by the single operator chain 
(SOC) (and successive improvements) approximation f20||2T|| . We have used the correlation 
corresponding to the Argonne v u+Urbana VII three-nucleon interaction model of Ref. [|2T] . 



A model with a weaker tensor force (the Urbana W14+TNI interactionl P^j ) was adopted for 
the CBF studies of Rl^ui) in Ref. ]TtJ and of the L/T spin responses in Ref. |f23 . 



The correlation given in Eq.(§) contains a q = 1, scalar (or Jastrow) component, f J (r), 
almost vanishing at short distances so that configurations where two nucleons are close 
enough are essentially suppressed in the wave function because of the short range NN re- 
pulsion. f J (r) heals to unity at large distances. The most important among the remaining 
components are, by far, the q = 5, spin-isospin, f UT (r), and the q = 6, tensor-isospin, f tT (r), 
ones, which are related to the one-pion exchange (OPE), long range part of the potential 
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The response is computed by considering up to correlated one particle-one hole (lplh) 
intermediate states. Admixtures of correlated particle (hole) states with 2plh(2hlp) ones are 



accounted for by the optical potential model of Ref. |24j| . This model has been proven to be 
fairly accurate in the momentum region scanned here, at least for the longitudinal response 
[TIJ. No correlated 2p2h intermediate states have been considered, since they contribute 



mainly at larger energies than those of the QE peak |17| , which is the focus of this work. A 



peak due the excitation of A-resonances is also present in the transverse response at high 
energies. The A peak is however well distinct from the quasielastic one and it has not been 
studied in the paper. It would have required the introduction of A-degrees of freedom in 
the nuclear wave function. 

The plan of the paper is as follows. The CBF approach to the transverse response 
is sketched and the current operators are briefly described in Section II . Section III is 
devoted to the discussion of the matrix elements of the current. The numerical results are 
presented and discussed in Section IV, together with a comparison with some of the available 
experimental data. Finally, some conclusions is drawn. 

II. THE CURRENT OPERATOR AND THE CBF THEORY OF THE 



TRANSVERSE RESPONSE 

The transverse response Rx(q,u!) is given by 

R T (q^) = W\m<l)\n}\ 2 6(u-u Jn ) , (4) 

A n 

where the sum goes over the intermediate excited states |n), having excitation energy 
u n , and j(q) is the electromagnetic current operator. j(q) has contributions from one-body, 
j^(q), and two-body exchange currents, j^(q). In the impulse approximation (IA) only 
j^(q) is retained. 

j^(q) is the sum of the convection, j^(q), and magnetization, j^(q), terms. However, 
it is known that the convection current contributes significantly to the non energy weighted 



sum of the response at very low q- values, where it becomes dominant p5[ . At momentum 
transfers > 150 MeV/c it provides only a few percent of the total sum, with a correction 
going as 1/q 2 . For this reason, (q) has not been considered here and we have approximated 



j (1) (q)~j$(q)=iG(?,")*r £ e 



e i=l,A 



i + r i>z l - n, 

V P o 



<*i x q , (5) 



where /j,q is the nuclear magneton and are the nucleon magnetic moments. The 
dipole parametrizion has been used for the nucleon form factor G(q,u), with a scale pa- 
rameter of 839 MeV. Studying the dependence of the response to the different available 
parametrizations of the electromagnetic nucleon form factors is beyond the scope of the 
present paper. 

jm (q) can be written in terms of the isoscalar and isovector transverse spin fluctuation 
operators, p^~ 0,1 (q), 

p°%) = i= £ e^fa X q) (r itZ ) , (6) 

V 1 i=l,A 

and, correspondingly, the IA electromagnetic response is related to the spin structure 
functions, Sy^q, u), by 

R I Aq,uj)=(-^=) IGMl'^lSl^^+^M^)} , (7) 



where p± = p p ±p n . Sy^.(q, u) have been computed in Ref. in FHNC/SOC. So, the 
details of the CBF calculation of the IA transverse response will not be discussed here. 
Most of the studies on exchange currents are based on simple meson exchange mechanims, 



in a non-relativistic description of the nucleon-meson-A systems |[26|| . In this case MEC 
associated with one-pion exchange are given by the sum of three terms, obtained by the 
corresponding Feynman diagrams: the contact, or seagull term, 3cont( c i)^ ^ ne pi ornc ) or pi° n 
in flight term, j^(q), and the A term, j^(q), associated with the excitation of intermediate 
A resonances. Analogous contributions arise from one-rho exchange. Their expressions can 
be found in several papers, and, among them, in Ref. 0]. These currents are consistent 
with NN potentials based only on one-boson exchange (OBE), so the continuity equation, 

V-j(x;r 1 ,r 2 )+z[u(r 1 ,r 2 ),p(x)] = , (8) 

with realistic potentials is not satisfied. 



In Ref. ||18|| , the Authors have derived exchange currents satisfying the continuity equa- 
tion with Argonne and Urbana Vn potentials. The current was separated in the sum of 
a model independent (MI) component, constructed from the NN interaction, and a model 
dependent (MD) one, purely transverse and given by the previously introduced tt and p 
MEC with intermediate A excitations. 

The most important MI currents are those associated with the isospin dependent parts 
of the NN interaction. The corresponding current operators, jps(q), jy (q) and jys(q), are 
given in Ref. ||18|| . They are expressed in terms of the Fourier tranforms of the v T {r), v aT (r) 
and v tT (r) parts of the potential |18j and reduce to the standard tt and p MEC if the OBE 



interaction is used. Argonne v\± has been used to generate the MI currents: they are close 
to those due to tt and p exchanges, even if the interaction does not strictly have the OBE 
form. 

The MD part of the current was associated to tt and p MEC with intermediate A ex- 
citations and to pTc-f and unry currents, all of them treated within the OBE model. The 
last two currents will not be considered here. Later, As were explicitly included in the wave 



function ||27|| , instead of introducing effective two-body operators acting on nucleon coordi- 
nates, as in standard first order perturbation theory. These components were generated by 
transition correlation operators, Uj. R = U^J A + U^ N + , acting on a realistic nuclear 



wave functions and obtained via the Argonne ^28 JT2[] potential, which contains A degrees of 
freedom. Hence, N— > A and A — > A transition operators were introduced in the one-body 
current. The A contribution to low-energy electroweak transitions was found to be smaller 
than that of perturbative theories as a consequence of this more realistic approach. It was 
also noticed that a first order perturbative treatment, consistent with that of Ref. |2^| , is 
obtained by using 



Upt (ra) = , 9 

m N - m A 



for the N— > A spin- and tensor-isospin transition correlations, where 

v a{t)rII are the 

corresponding transition components of the Argone t>28 potential. 

The intermediate A current generated in this way will be used. Its configuration space 
expression 



y i<3 
.trlli^t. ^ CM (~ w - \ \(„.crrll( m \ n ,trII 



+ 3upT (r)r x q(tr i • f)] - (n x tj) z [(u P T ^ (r) - u% (r))((Ji x aj) x q 
+ 3ii^iXr)xq(a r f)]} + i#i] . (10) 



r = Tij and /i 7 7VA = 3/x is the transition magnetic moment adopted in Ref. p7| , whose 



value is ~ 30% smaller than the static quark model prediction p9| . For the A form factor 
we take 

G.(^)=(l + |)"(l + |)" /2 , (ID 
where the cutoff masses are A x = 1196 MeV and A 2 = 843 MeV. 

As a result of the separation of the current into one- and two-body pieces, the response 
is given by 

R T (q, u) = R'/iq, oj) + R^ n \q } u) + R^ IEC (q } u) , (12) 



where Bjp lt (q,u) is the interference contribution between j^^q) and j^(q), 

i?^(g,c) = i5:[(0|j( 1 )(q)|n)(n|(j( 2 )(q))t|0) + c.c.]^-c n ) , (13) 
A n 

and 

< £C (g^) = 7El(o|j (2) (q)Hr^-^.) • (14) 

A n 

Rj! lt (q,uj) and the leading term in R^ IEC \q, uj) only have been computed. 

Let's now shortly discuss the most important aspects of the CBF perturbative expansion. 
The goal of using a correlated basis, embodying directly into the states some of the relevant 
physical effetcs (as the short range repulsion), is to obtain a rapidly converging expansion. 
The obvious price to be paid lies in the greater difficulty in evaluating matrix elements, 
even at the zeroth order. However, nuclear matter studies demonstrated that CBF based 
perturbation theory is actually fast converging, in the sense i) that for many observables 
the zeroth order is already a good estimate and ii) that the inclusion of the first CBF 
intermediate states is often sufficient to provide a quantitative agreement with the empirical 



values. This happens for ground state properties, as the energy |30| and the momentum 



distribution ||31|| , as well as for the already mentioned longitudinal and spin responses and 



the nucleon spectral function B3|. The inclusive electron scattering cross section at high 



momentum transfers has been computed in Ref. ||33|| , using the CBF spectral function, and 
satisfactorily compared with that extrapolated from data on laboratory nuclei. 

Hence, relying on these facts, we will compute Rt{q,uj) by inserting in the intermediate 
states summation of Eq.([|) only correlated lplh states 



|ph) = S 



i<j 



\ph.) FG . (15) 



This choice is also justified on the ground that we are interested mainly in the QE peak. 
The quantitative study of the large energy region would have required the insertion of higher 
excited states with both nucleonic (as 2p2h) and explicit A isobar excitations (as A-h). 

The lplh transverse response is then given by 
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R¥ l \q,u) = jY<\m<l)\ph)\ 2 5(uj-e p + e h ) , (16) 

ph 

where e x =p,h — ( X I-^1 X ) ~~ (0(^Z"|0) is the CBF variational (or zeroth order) single particle 
energy of the x-state. 

The \p\h response has sharp energy boundaries, qualitatively like the Fermi gas, ruled 
by the real part of the CBF optical potential, Uk = &k — h 2 k 2 jlm^ f3(|. The decay of Iplh 



states into 2p2h ones has the main effect of introducing a large energy tail and redistributing 
the strength. In the perturbative expansion this is accounted for by self-energy insertions 
on top of the particle or hole line. The real part of the self-energy provides a perturbative 
correction, Se x , to the variational single particle energy and the imaginary part induces the 
spreading of the response to high cu-values. The microscopically computed CBF self-energy 
was used in the longitudinal response calculation of Ref . |10| . In Ref . |24| these corrections 
were estimated by folding the lplh response with a width W(u) given by the imaginary 
part Wq(u) of the optical potential divided by the nucleon effective mass, 

with Wq(u) ~ llc<j 2 /(4900-|-ct; 2 ), in MeV. This procedure, numerically much less involved 
and equivalent to retain the on-energy shell part of the self-energy only, was checked to be 



reliable in the momentum region of interest [|I0] and, for this reason, has been adopted here. 



III. CBF MATRIX ELEMENTS OF THE CURRENT 

This Section will focus on the main features of the MEC CBF matrix elements. As far 

• (2) 

as the MI currents is concerned, only the leading jps(q) will be discussed. However, the 
contributions of the less influent (q) and jy^(q) have been evaluated and will be presented 
in the next Section. 

(2) 

The configuration space PS current is written as a sum of two terms, ipsci.^) an d 
jp i 5 7r (q). The first one coincides with jclnti*!) f° r the one-pion exchange potential, the 
latter with j^ 2 ^(q). Their expressions are 

JP5,c(q) = C PS E 3(r* x TjUe^gpsir)^ ■ r) + i # j} , (18) 

i<j 
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and 



jfs,M = CpsT.Hn x ^,e^^|^Ma, • r) + a fa • f) + fa • a,)] 

%<j 

, GpsM Gps, 3 (r) , s 
+ « CilCj • q) - i vfai ■ q) 

- i fai ■ r)(a\j ■ q) + 2 r(or 3 - • r)(a< ■ q) + 

- G PS)6 (v)v (a, ■ qfa ■ q) - «^^f (a, ■ r)(a, ■ r)} . (19) 



The functions gpsij') and Gps, a =i,7( r ) are defined in Ref. fl8fl , R = + ij)/2 and 
Cps = G(q,u). 

In CBF the non diagonal matrix elements, £ps,c*(ti-) (q; P, h) = (Oljp^^q^p, h) are 
computed by cluster expansions in Mayer like diagrams, built up by dynamical correlations, 
y(<j)y(p) _ S ql S p i, and by exchange links. Infinite classes of cluster terms containing Jas- 
trow correlations are summed by FHNC. Less massive summations, similar to SOC, can be 
performed for the operatorial correlations [I0"|,21|. 



In the case of simple Jastrow correlated wave functions (f q>1 = 0), (,p S ,c * s gi ven by 

^P5,c(q; p ; h ) = Sq-p+tfCps 



D(p)D(h) 

Po J drgp S (r)g cc (r) {e* p - r (r - lo x r) + e^ r {r + taxv)} , (20) 



where t z = r Pj2 = Th, z , & = & P = &h, D(x = p,h) = 1 — X cc (x) (see Ref. |24| for 
the definition of X cc ), g cc (r) is the exchange FHNC partial radial distribution function 
fl9|| and po is the nuclear matter density. At the lowest order of the cluster expansion, 
D(x) = 1 and gcc{ r ) — ~ [/ J ( r )] 2 ^(^P r )/4, where kp is the Fermi momentum and l(x) = 
3[sin(x) — xcos(x)]/x 3 is the Slater exchange function. So, Jastrow correlations suppress the 
short range part of the pion-like exchange interaction. This effect is visualized in Fig.|l], where 
the Argonne v 14 gps and the FHNC g cc functions are shown. Moreover, the Figure compares 
g cc (FHNC) with its uncorrelated Fermi gas counterpart, g cc (FG) = — Z(fcpr)/4, to stress 
their different short distance behaviors. Finally, the simple 7r exchange gps(Tr)(f%/4ir = 
0.081) is given. As it was already stated, it results to be close to the full gp$- 
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Eq. pO|) sums all cluster diagrams factorizable in products of dressed two-body diagrams. 
Non factorizable diagrams involving three particles have also been taken into account, even 
if they do not appear in the expression. 

Operatorial correlations were introduced in Ref. |23[] for the spin responses within both 



the dressed two-body (D2B) and the SOC approximations. In the first approximation cluster 
diagrams containing Jastrow correlations are summed to all orders, while contributions 
from the other components are evaluated at the two-body level (it corresponds to the Wo 
approximation of Ref. |20fl). D2B came out to be very accurate and it has been used here 



to evaluate the IA response. A linearized version of this approximation (D2B/L), including 
only contributions linear in f (JT,tT , has been adopted for the MEC matrix elements. However, 
in some cases quadratic terms, as well as those containing the other operatorial components, 
have been computed and checked to give small corrections. 

Following the above scheme for the treatment of the q > 1 correlations, £ps,c has an 
exchange part, obtained by the replacements 

gcc(r)i - 9cc {r)v jl + 2 rr( ^7 (r 4 ) fT(r) I , (21) 

f f^ir) - f tT (r) 1 
g cc {r)t{o x r) -> 9cc (r)i(a x r) |l + 2 J - K J f J K } \ , (22) 

in Eq . (p0|) , and a direct one 

£fs,c(q; p. h ) = <* q -P+h*Cps 



D(p)D(h) 

r r (r)-f tr (r) 



p J drg PS (r)[g dd (r) + g de {r)\e^2i{a x f) | J [r ' | , (23) 

where gdd,de( r ) are the direct-direct and direct-exchange FHNC partial distribution func- 



tions. 

;(2) 



jpsc(q) is the leading two-body current and it is the only one used to estimate 
B^ EC {q,u) of Eq.H 
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jpg 7r (q) is the sum of several components. Following an obvious notation, we will refer 
to the corresponding parts of its matrix element as 0^s>(G\ps>=i,7)- They are given in the 
Appendix. 

In presenting the matrix elements of the jg D (q) current, £ MO (q; P, h) = (0|jg ) D (q)|p, h), 
we introduce, for the sake of brevity, the notations 

C A = lG A (q,cj)ii jNA , tf(r)=vgf I {r)-v& I {r) , f 2 A (r) = 3«#'(r) . (24) 

£,md is given by 



£ MD (q; P, h) = £gb(q; p, h) + Cmd(^ P, h) 



exc/i / 



(25) 



where 



Q T d(^ P> h ) = ^- p+h iC A q 



4r, 



Po / dr[g dd {r) +^ e (r)]e iqr 



+ / 2 A (0 



cr X q) ( 1 
cr x q 



D{p)D{h) 

V aT {r)+f tT {r)' 



f J (r) 

2f r W-/* T (r)' 



— 3(cr ■ r)(f x q) 



/ J (r) 



+ (cr • r)(r x q) 1 + 



/ tT (r) 
f J {r) 
f tT (r)-r T (rY 
f J {r) , 



4r 2 f 

8q- P +hiC A q I =po / dr[g dd (r) + ^ e (r)] [(cr x q) + (a ■ r)(r x q) 

. / n 0n\ r>(h\ J 



D(p)D(h) 



f J (r) 



f J (r) 



(26) 



Gd (q; p> h) = 5 q _p +h 2C A g 



2tv 



D(p)D(h) 



p / drflr^r 
— 9(a • f)(r x q) 



e tp.r + e «h-r 

/ tr (r) - 



/ A (r) 



V x q) 4 + 3 



/ J (r) 



+ / 2 A M 



(ffxq) 1 + 2 



/J(r ) 

r r (r) - f r (r) 
/ J (r) 



+ (<r-r)(r xq) 1+6 



3r r (r)-/^(r) ' 
/ J (r) 



(27) 



It is worth noticing that the second integral in the r.h.s. of Eq.(p6|) gives a correction 
independent on the modulus of q and vanishing for uncorrelated wave functions as well as 
for Jastrow correlated ones. 



13 



Fig.0 compares the f^ 2 functions obtained with the Argonne f 28 potential, f^ 2 (A 28 ), 
with those from the OPE potential, f^ 2 (^)- They differ from each other at short distances, 
and are therefore expected to provide similar results in a correlated model. 

IV. RESULTS AND DISCUSSION 

In this section we present and discuss the inclusive transverse electromagnetic response 
of symmetric nuclear matter at saturation density, po — -16/m~ 3 or kp = 1.33/m -1 , as 
obtained in CBF theory. We recall that the correlation operator contains Jastrow, spin, 
tensor and isospin components. The correlations are fixed by minimizing the ground state 
energy, computed in FHNC/SOC and using the Argonne t>i4+Urbana VII three nucleon 



interaction ||34|| . The SPR model for the two-body current operators of Ref.s flqj27|1 has 
been used and the one-body current is given by the magnetization part alone, neglecting the 
convection term. The response is computed including correlated lplh intermediate states 
and then folding Rj? lh (q,Lj) with a parametrization of the on shell imaginary part of the 
CBF optical potential. 

We begin by studying the effect of the Jastrow, short range correlations and comparing 
with the free Fermi gas. The lplh responses in both models are given in Fig.^ at q — 300 
(a), 400 (b) and 550 (c) MeV/c. In the Figures the one-body response (IA), the interference 
terms between one-body and jpsci^) (OB/C), one-body and jp5. 7r (q) (OB/V), one-body 

(2) (2) 

and j MD (q) (OB/A) and the quadratic jpgc?(q) (C/C) term are shown, together with the 
resulting total response (TOT). 

The shift of the Jastrow responses to higher energies respect to FG is due to the use 
of the CBF real part of the optical potential; moreover, as in the case of the longitudinal 



response [[UJ , Jastrow correlations quench the quasielastic peaks in absolute magnitude. The 



computed quadratic term is always negligible and will be disregarded in the remaining of 
the paper. Two interference responses (OB/C and OB/ A) have been evaluated also for the 
one pion exchange currents. The results at q — 400 MeV/c are displayed in Fig.|3|(b) as up 
and down triangles, respectively. They are almost coincident with the SPR responses, as the 
differences in the short range behaviors of the A currents are washed out by the correlation. 
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The numerical accuracy of the FG responses has been checked against the analytical 



approach of Ref. |T7| finding complete agreement. So, an overall reduction of the lplh QE 
peak due to MEC's is confirmed in the FG model and found even in the Jastrow correlated 
case. 

As we have already mentioned in the Introduction, this result is in sharp contrast with 



realistic estimates of the transverse response in light nuclei |yj and also with some recent 



RPA calculations in 12 C [35], both of them pointing to a ~ 20 — 30% increase of the strength 



in the quasielastic region due to MEC effects. We stress once again that this conflict is not 
resolved by the inclusion of state-independent, scalar, short range correlations. 

In Ref. |23| it was shown that the transverse spin responses, making up R^(q,u), may 



be affected by the non- Jastrow correlations. In particular, the biggest effect was found in 
the isovector S^ 1 , as the leading correction is proportional to the large tensor- isospin corre- 
lation, f tT (r). Hence, the correlation operator (|3]) has been used to estimate the corrections 
to the Jastrow response due to the non scalar components, in the D2B/L approximation. 

Fig.[|(a) reports the correlation operator results for the interference lplh responses at 
q = 400 MeV/c. We find that, at the QE peak, the operatorial correlations quench the 
OB/C and OB/V responses respect to the Jastrow case. The use of the D2B in place of the 
D2B/L approximation does not change appreciably the outcome. The quenching is more 
pronounced for the OB/7T term, where a positive tail is added at large energies. 

The effect is dramatic in OB/A, as the correlation operator correction largely cancels 
the Jastrow response, yielding a positive net result. The origin of this cancellation is found 
in the tensor-isospin correlation contribution to the second integral in the r.h.s. of Eq. (p6|) . 
In fact, the OB/A response obtained by setting f tT = in the integral is much closer to 
the Jastrow curve, as it is shown in the Figure by the x signs. The convergence of the 
cluster expansion has been checked by computing the integral i) in D2B approximation and 
ii) adding Jastrow dressed three-body non factorizable diagrams, linear in the operatorial 
components of the correlation. The result is practically indistinguishable from the D2B/L 
response and it is not given in the Figure. 
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A similar (and even more enhanced) behavior is found at larger momenta. In Fig.[|(b) 
we show the interference terms at q — 550 MeV/c. Now the OB/C response does not 
significantly differ from the Jastrow estimate, whereas OB/V and OB/ A are largely modified 
by the insertion of the operatorial correlations. In particular, we obtain a large, positive 
OB/ A component in place of a large, negative Jastrow counterpart. At this momentum 
transfer and using only long range RPA correlations, the Gent group seems to find a more 
negative OB/7T response and an almost vanishing OB/ A one |36 . 



The sensitivity of the OB / A response to the shape of f tT has been examined by using 
the correlation operator corresponding to the Urbana W14+TNI model. This interaction has 
a weaker tensor component, producing a smaller tensor correlation. The results for OB/ A 
are displayed in Fig.s| as full circles and point to a clear dependence on the details of this 
correlation. 

To conclude the analysis of the lplh response, we give in Figures |5] the total Rj^ lh (q, u) 
together with the IA estimate. The total (TOT) response includes interference contributions 
also from the j[/*(q) and jy^(q) currents, whereas the curves labelled as TOTps only contain 
jps(q) (as the Jastrow, TOTj, points). MEC provide an enhancement of the IA response 
in the QE region ranging from ~ 20% at q = 300 MeV/c to ~ 10% at q = 550 MeV/c. The 
enhancement is due to the presence of tensor- and isospin-dependent correlations; scalar 
Jastrow correlations actually produce MEC contributions quenching the corresponding IA 
values. The exchange of p-like vector mesons, originating the V and VS currents, results in 
an additional enhancement of Rj? lh (q,uj), which, however, decreases with q. 

Finally, we compare the transverse NM response with some of the available experimental 
data from real life nuclei. Rr(q,uj) is obtained by folding the computed Rj? lh (q,Lj) with 
the imaginary part of the optical potential, as described in Section II. Figures ^| compare 
A x R T (q,uj) of NM in IA and IA+MEC with data from Ref. § (x), Ref. § (circles) and 
Ref. (black circles) for 40 Ca. The Bates data (Ref. §) actually refer to q = 330 MeV/c 
in Fig .0(a) and q = 400 MeV/c in Fig.^(b), as the Saclay data (Ref. ||) in the same Figure. 

The discrepancies between the Saclay and the Bates and Jourdan (Ref. |§) results are 
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evident and already stressed in previous literature. It must be noticed that the last two 
sets of data appear to be compatible among each other. This fact gives confidence in the 
Jourdan's L/T separation procedure, which employs world data on inclusive quasi-elastic 
electron scattering as obtained by different experiments. 

The NM responses are closer to the results of Ref.s |||9|] than to those of Ref. Q and the 
experimental QE peak lies between the IA and the full responses. At q = 300 MeV/c IA 
seems to better reproduce the peak, whereas, when moving to higher momenta, the inclusion 
of MEC improves the agreement with the experiments. 

A similar trend is found in 56 Fe, as shown in Figures |7]. 

V. CONCLUSIONS 

The inclusive electromagnetic transverse response of symmetric nuclear matter has been 
evaluated within the correlated basis function perturbation theory. The adopted correlation 
operator has a scalar, short range component and important tensor- and isospin-dependent 
parts, its structure being similar to that of realistic nucleon-nucleon potentials. Both one- 
body and two-body meson exchange currents are considered, the latter in a model which 
satisfies the continuity equation with the Argonne vu potential and contains intermediate 
A-isobar excitation currents. Ground and lplh correlated states are included and the decay 
into 2p2h states is implemented by folding B)f lh (q, uo) with the imaginary part of the optical 
potential. 

Our results indicate that MEC, evaluated in a Jastrow correlated model, quench the IA 
response. In this case, the situation is qualitatively close to what was found by Amaro et 



al. in Ref. [17]] in both shell and Fermi gas models. The net quenching mainly originates 
from a strong cancellation between the positive contact and the negative A terms. 

The introduction of tensor-isospin-dependent correlations drastically changes this pic- 
ture. The A contribution is largely modified, as it becomes positive and increasing with the 
momentum transfer. As a result, MEC produce an extra-strength (10-20%) in the QE peak 
region. This is in agreement with exact GFMC calculations in light nuclei. 

p-like exchange currents give a small additional enhancement. We also found that using 
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standard one-boson exchange currents does not significantly change our results. 

Two recently derived experimental responses in 40 Ca have consistently lowered the QE 
peak respect to previous estimates. The new data and the CBF NM responses are in 
reasonable agreement and the comparison seems to show too large MEC effects at low 
momenta. The obvious caveat to bear in mind is that this comparison is made between 
finite nuclear systems and infinite, homogenous nuclear matter. CBF has been recently 
extended |37] to deal with ground state properties of nuclei as heavy as 208 Pb, with Jastrow 
and isospin-dependent correlations. It is conceivable that, in the near future, will be possible 
to use the theory to microscopically compute the finite nuclei responses, employing richer 
correlations, as those of nuclear matter. Presently, the density dependent NM results might 
be used in local density approximation for a closer comparison with the experiments. 



Moreover, relativistic corrections could affect the response, both for the IA and MEC j38 
as well as the inclusion of explicit A degrees of freedom in the nuclear wave function. This 
last topic may be addressed within CBF theory, particularly if one wants to quantitatively 
study the dip and A peak regions. In this respect, it will also be necessary to consider the 
contribution from 2p2h intermediate correlated states. Work along these lines is in progress. 
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APPENDIX 

In this Appendix the particle-hole non diagonal matrix elements of the components of 
jp^ 7r (q) are given in D2B/L approximation: 



£psJn(GpS,l) — ^q-p+h^CpS 1 



D(p)D(h) 
poj dr^^-[g dd (r) + g de (r)]e^ 2 

6t, 



6* * x v)Ufl 
f J (r) 



(28) 



D(p)D(h) 



f , GpS,l( r ) l \ «Q-r 

po J dr g C c{r)e " 



f J (r) 



f J (r) 



6t z 



Po J dr 



D(p)D(h) 
G PS ,2(r) + Gp5 i3 (r) 



[g dd (r) + g de (r)]e^ 2 



f tT (r) f aT (r) 
[3{q-r){axr)-3i(axr)-q-2(axq)]^-jy--2{axq) J 1 ; 



f J (r) 



f J (r) 



Cp X sl(GpS,2+z) — 5q- p+ hlCps 



6tv 



^D{p)D{h) 



Po J dr 



G PS ,2( r ) + G PsA r ) 



gcc{r)e 



iQr 



(cr x q) - 3[(q • f)(cr x r) + (cr ■ r) (q x r)] 



f T (r) 
/ J (r) 



6r 2 



po y 



D(p)D(h) 
Gp SA (r) + G PS , 5 (r) 



[to(r)+^ e (r)]e^ r / 2 



f (cr x f ) • q 



f J (r) 



Cp x s h 7r(G ps, 4+5) - Sci- P +hiCps 



6r 2 



D{p)D{h) 



Po J dr 



GpsA r ) + G PS , 5 (r) 



g C c(r)e lQ r [r(a x r) ■ q] , 



^P5,7r( G P5,6) - , 



Cp x sAGps,e) - -5 q - p+ h?Cp5 



6tv 



^D{p)D{h) 



po J drG PSfi {r)g cc {r)e^ r 



3«r(q • r)(cr x f) • q 



f r (r) 



Kps,AGps,i) = t,ps,n\.Gps,7) = U , 



exch , 



where Q = (p + h)/2. 
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FIGURES 

FIG. 1. gps functions for the Argonne vu, gps(Ai4), and the ir exchange, gps(n), potentials. 
g cc (FHNC) and g cc (FG) are the correlated and uncorrected exchange functions, respectively (see 
text). 

FIG. 2. /f^ functions for the Ar gonne V2s, /i^O^s)) and the tt exchange, f^i^)^ potentials. 

FIG. 3. lplh transverse response at q = 300 (a), 400 (b) and 550 (c) MeV/c for the FG and 
Jastrow models. The Figure shows the partial (one-body, interference and one quadratic) and 
the total responses. The interference terms for the contact (triangles up) and A (triangles down) 
currents in the one pion exchange model at q = 400MeV/c are given. 

FIG. 4. lplh interference responses at q = 400 (a) and 550 (b) MeV/c in the correlation 
operator model and comparison with the Jastrow model. See text. 

FIG. 5. Total lplh responses at q = 300 (a), 400 (b) and 550 (c) MeV/c in the correlation 
operator model and comparison with the Jastrow model and the IA responses. See text. 

FIG. 6. Transverse responses at q = 300 (a), 380 (b) and 570 (c) MeV/c for 40 Ca and nuclear 
matter. See text. 

FIG. 7. Transverse responses at q = 380 (a) and 570 (b) MeV/c for 56 Fe and nuclear matter. 
See text. 
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